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Abstract 

We study open B-model representing D-branes on 2-cycles of local Calabi-Yau geometries. 
To this end we work out a reduction technique linking D-branes partition functions and multi- 
matrix models in the case of conifold geometries so that the matrix potential is related to the 
complex moduli of the conifold. We study the geometric engineering of the multi-matrix models 
and focus on two-matrix models with bilinear couplings. We show how to solve this models 
in an exact way, without resorting to the customary saddle point/large N approximation. The 
method consists of solving the quantum equations of motion and using the flow equations of the 
underlying integrable hierarchy to derive explicit expressions for correlators. Finally we show how 
to incorporate in this formalism the description of several group of D-branes wrapped around 
different cycles. 
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1 Introduction 



Singular Calabi-Yau spaces are more and more frequently met in string compactifications. The 
reason is mostly the fact that compactifications on regular Calabi-Yau spaces do not seem to be 
able to describe crucial features of realistic physics. On the contrary the presence of a conifold 
point, PP, in a Calabi-Yau opens new prospects: in conjunction with fluxes and branes it may 
allow for warped compactifications, which in turn may create the conditions for large hierarchies 
of physical scales. On the other hand singular Calabi-Yau's with conifold singularities seem to be 
necessary in order to realize low energy theory models with realistic cosmological features. The 
hallmark of a conifold is the possibility of resolving the conifold point in two different ways, by 
a 2-sphere (resolution) or a 3-sphere (deformation). This leads, from a physical point of view, 
to a geometric transition that establishes a duality relation between the theories defined by the 
two nonsingular geometries (gauge-gravity or open-closed string duality), |2| Ej. In summary 
conifold singularities are at the crossroads of many interesting recent developments in string 
theory. It has therefore become customary to study theories defined on conifolds, i.e. singular 
non compact Calabi-Yau threefolds, as calculable and well-defined models to approximate more 
realistic situations. 

Given the crucial role they play it is of upmost importance to find methods of calculation 
for theories defined on conifold geometries. In this sense two main tools have been devised: 
topological field theories and matrix models. Topological field theories are truncations of full 
theories: one gives up the knowledge of the dynamical sectors of a given theory, drastically 
simplifies it by limiting it to the topological sector and ends up with a theory where very often 
many quantities (correlators) can be explicitly calculated. However even topological field theories 
are sometimes not easily accessible to explicit calculations. Here come matrix models to the 
rescue. Sometimes, like in the examples of this paper, topological field theories can be shown to 
be equivalent to matrix models. This makes life easier, especially when the matrix models have 
couplings of special type. In this case one can rely on the integrable structure underlying the 
model (the Toda lattice hierarchy, [55]) which usually provides algorithmic methods to obtain 
the desired results in an exact and controlled framework. The case of matrix models with more 
general couplings is more complex and represents a challenge that people have started to tackle 
only very recently. 

In this paper we would like to elaborate on an idea that has recently received increasing atten- 
tion: how data about the geometry of a local Calabi-Yau can be encoded, via a topological field 
theory, in a (multi-) matrix model and how they can be efficiently calculated. The framework we 
consider is IIB string theory with spacetime filling D5-branes wrapped around two-dimensional 
cycles. This geometry defines a 4D gauge theory, (111310 El- On the other hand we can consider 
the open topological B model representing the strings on the conifold. The latter has been shown 
by Witten long ago to be represented by a six-dimensional holomorphic Chern-Simons theory, 
1^1. When reduced to the two-dimensional cycle this theory can be shown to boil down to a ma- 
trix model. In particular, if we wish to represent the most general deformations of the complex 
structure satisfying the Calabi-Yau condition, we end up with very general multi-matrix models. 
This point of view was advocated in ,7,. In this paper we concentrate on the topological string 
theory part of the story, and ignore both the 4D gauge theory part and the closed string theory 
side, which is attained by shrinking the resolved sphere to a point and passing to the deformed 
picture in which the singularity is replaced by a three-sphere (as would be possible at least for 
the cases corresponding to one and two-matrix models). 

As already pointed out the general idea underlying our paper has already been developed in 
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a number of papers, jlUl 1131 ITTl 1141 112j . Here we would like to concentrate on particular aspects 
that have not been stressed or have been left aside in the previous literature. The first question 
we concentrate on is the reduction from the six-dimensional holomorphic Chern-Simons theory to 
a two-dimensional field theory. We wish to understand what degree of arbitrariness this passage 
implies, so as to be able to assess whether the information we gather from the reduced theory 
is intrinsic or depends on the reduction process. Our conclusion is that the reduced theory does 
not depend on the reduction procedure. 

The second point we deal with is whether there are limitations on the general form of the 
potential we find for the multi-matrix model. We do find some conditions although rather mild 
ones. Finally we concentrate on the subclass of matrix models represented by two-matrix models 
with bilinear coupling. In this case the functional integral can be explicitly carried out with 
the method of orthogonal polynomials. We show, using old results, how one can find explicit 
solutions: this is done by solving the quantum equations of motion and utilizing the recursiveness 
guaranteed by integrability. All the data turn out to be encoded in a Riemann surface (plane 
curve), which we call quantum Riemann surface in order to distinguish it from the Riemann 
surface of the standard saddle point approach. In particular we are able to prove that the exact 
solutions found in this way are more in number than the ones found by the saddle point method. 

The paper is organized as follows. In the next section we calculate the reduction of the 
B-model open string field theory corresponding to the wrapped D-branes on the 2-cycle of the 
conifolds and how it depends on its complex moduli. The result is given in terms of a family of 
multi-matrix models. In section 3 we discuss explicitly what multi- matrix models we do get and 
we draw the geometrical engineering scheme for their realization with D-branes. In particular 
we show how to obtain two-matrix models with bilinear couplings. In section 4 we review some 
general properties of the above class of 2-matrix models, their integrability and the relative genus 
expansion. The section is concerned with particular method of solving them which is based on the 
quantum equations of motion and the flow equations. Section 5 contains several explicit examples 
of models solved with this method. We compare these solutions with the ones obtained with the 
usual saddle point/large N expansion. All these solutions are interpreted as describing the physics 
of A'^ D-branes wrapped around a 2-cycle. In section 6 we show how the physics of several group 
of D-branes wrapped around different cycles can be incorporated in the exact scheme proposed in 
this paper. Finally, section 7 is devoted to some conclusions and open questions. The Appendix 
extends the approach of Section 2 to local CY geometries around 4 cycles. 

2 Reduction to the branes and multi- matrix models 

In this section we show in detail how the reduction of the topological open string field theory 
(B model) to a 2-cycle in a local CY geometry is equivalent to a multi-matrix model whose 
potential is parameterized by certain deformations of the complex structure of the non compact 
CY space. Actually we will elaborate a more general framework. We consider three dimensional 
Calabi-Yau geometries built around a generic Riemann surface S of any genus which is then the 
non trivial 2-cycle we wrap the D-branes around. The normal bundle is specified by assigning a 
rank two holomorphic vector bundle V over T, and the CY condition constrains the determinant 
line bundle to equal the canonical line bundle of S, while V is otherwise generic. 

It is in this generic setup that we study the problem of the reduction of the holomorphic 
Chern-Simons action functional to the D-brane world volume. This depends on the (0, l)-part of 
a connection on a U{N) gauge bundle E which we take to be trivial. Due to the non triviality of 
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the geometry of the normal bundle, in order to specify the reduction mechanism, we will need to 
choose a trivialization of the bundle by a reference non-degenerate bilinear structure K over it. 
We will show that choosing the reference connection to be the (generalized) Chern connection 
of the bilinear structure K makes the overall result actually independent both on the reference 
bilinear form K and on the base representative of the 2-cycle in the total CY space. The resulting 
reduced theory is a generalized holomorphic 6-c (/3-7) system on S where the two bosonic fields 
span a section of V. These are minimally coupled to the reduced gauge connection. 

After the above preliminary construction, we consider the effects of varying the complex 
structure of the total space. Actually we will study constrained variations leaving the complex 
structure on E fixed and preserving the CY condition. These can be seen to be parametrized by 
a set of geometric potential functions on the double intersections which are the Cech cohomology 
representatives of the complex moduli we are varying. These, in the spirit of Kodaira and Spencer, 
can be used to parametrize local singular coordinate changes which specify the variation of the 
complex structure. 

In order to be able to explicitly deal with the moduli space of the conifold complex structures, 
we then limit ourselves to the genus zero case, i.e. S = P^. In this case, the cycle has a single 
complex structure and so the above analysis is enough to cover the full moduli space. Since 
the transverse fields are sections of the normal bundle V, the singular coordinate transformation 
defines the deformation of the reduced theory action in a well defined way which is parametrized 
by the geometric potential at the intersection of the north and south pole. We calculate this 
explicitly in the generic case of CY deformations of the 0{n) © 0{—n — 2) reference complex 
structure on the conifold. The result we find is that the partition function generically reduces 
to an (n+l)-matrix model whose potential is obtained from the geometric potential in a specific 
way. 

Let us stress that the above result with S = P^ and n = was presented in [3] where it was 
suggested that it can be obtained by refining a sketchy calculation in JS]- The result presented 
here is a generalization thereof for Riemann surfaces of arbitrary genus and, for arbitrary n, on 



The same reduction method can be applied also to non compact local CY geometries build 
around a 4-cycle. In such a case the CY condition specifies the normal line bundle to be the 
canonical line bundle on the base complex manifold. In the appendix, we elaborate it in a generic 
case and find the reduced holomorphic Chern-Simons theory in the form of an holomorphic BF 
model (see [16j for recent discussions on this model). It was recently suggested in [TH- that it 
describes the topological open strings for D-branes wrapped around the above four cycles. 

2.1 Reduction in the linear case 

Let us consider non compact six dimensional geometries built around a Riemann surface S as 
the total space of a rank 2 holomorphic vector bundle V with GL{2, C) structure group. 

Any atlas {^/(q,)} on S extends to an atlas on CY{Ti, V) by U(^a) = f^(a) >^ C^. The complex 
manifold is defined by the overlapping conditions 






in any double patch intersection ?7(^) n C/(^). 
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Requiring the complex manifold to be of the Calabi-Yau type, restricts det V to be equal to 
the canonical line bundle T^^''^\T,) so that under this condition the total space of V is equipped 
with the holomorphic (3, 0)-form i7 = A dw^ A dw"^, where z is a local coordinate system on S 
and on the fibers. This condition is just detMj-^-j^^) x f(^a){^) ~ ^ ^^'^ ^® consistent with 
triple intersection conditions. We denote this manifold CY{T,,V). 

Let us consider the topological open B-model on CY{Ti, V), which can be obtained starting 
from open string field theory, Because of the drastic reduction of the degrees of freedom 
due to the huge gauge symmetry present in the topological string, the string field theory is the 
holomorphic Chern-Simons (hCS) theory on Cy(S,y) for a (O,l)-form connection on a U{N) 
bundle E, where N is the number of D-branes wrapped around S. For simplicity, we will restrict 
to the case in which E is trivial. The action of hCS is 

S{A) = —[ C, £ = ATr ( ^^Aa^+ ^^A^A^ ) (2.2) 

9s JcY{Y.y) V2 3 J 

where A G T^^'^) (Cy(S, V)). 

The dynamics of D-branes wrapped around the 2-cycle S can be described by reducing the 
open string field theory from the total space manifold to the D-brane world-volume. Since 
the bundle V is non trivial, in order to properly define the brane theory, the reduction of the 
lagrangian has to be coherently prescribed patch by patch by a trivialization procedure in such 
a way that the end product is independent upon the particular trivialization we use. 

As it is evident, the embedding equations for T, in CY{T,, V) are just = 0. In any local 
chart, the fibering structure defines a local notion of parallel and transverse directions along 
which we split A = Azdz + AidiB'^. The parallel part Az glues on double patches intersections as 
an invariant (O,l)-form on S only when restricted to the base, while otherwise gets also a linear 
contribution in w due to the generic non triviality of V The transverse coefficients Aj glue* as a 
section of V* . 

Because of this, since the reduction to the base has to be performed covariantly, we have 
to expand As = - A-^T^--.w~^ and Ai = Aj, where A^dz E r(°'i)(S), A- G V* and dzT^--. 
is the (0, 1) component of a reference connection of V. The reduction process is defined by 
specifying the subfamily of A connections we limit our consideration to. Our prescription is that 
the matrix valued dynamical fields {Az,Ai) that survive the reduction are those independent 
of the coordinates along C^. A direct calculation from the lagrangian C in ()2.2() for the above 
reduced configurations gives 

L = n A Tr Q [AjDzA-j + ^^rl^^}^ du/ Adz A dvP (2.3) 

where is the covariant derivative w.r.t. the gauge structure. Notice that the above does not 
depend on the base representative, that is on the values of . 

Another way of justifying the reduced lagrangian 1)2. 3|) is the following. We start from ()2.2|) 
and replace the exterior differential d by the covariant differential D = 5 -|- F on Cy(E, V) and 
impose that "^jAj = = DjAz- The latter conditions are satisfied as follows. Due to the local 
product structure of CY(T,, V), we can suppose without loss of generality that the only nontrivial 
component of the connection relevant to the problem are F^^ = F^_. Therefore TijAj = simply 

ZJ JZ J 

means that Aj does not depend on w^, while TijAz = can be integrated and leads precisely the 
expression for Az given above. 



"We denote by V* the dual vector bundle, glueing with {M and by V the complex conjugate one. 
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Ea. ()2.3() represents a six-form. Our purpose is to restrict it to a two form defined on S. This 
has to be consistently prescribed patch by patch. That is, we have to couple the above reduction 
with the contraction of the differentials along the fiber directions to obtain a well defined (1, 1)- 
form on S. To define this operation, let us consider a bilinear structure K in V , that is a local 
section K £ T{V <^ V), the components K'^^ being an invertible complex matrix at any point. 

The derivation of the basic (1, l)-form is realized patch by patch with the help of K as con- 
traction of the hCS (3,3)-form Lagrangian by the two bi- vector fields k = ^CijK^^K^'^-^-^p^ and 
p = 5^*"' gf^Tgfjj- Notice that k S detV and p £ detV* = {detV)~^ so that the combined applica- 
tion of the two is a globally well-defined operation. Calculating then the pullback Lagrangian, 
we obtain 

Cred = ipAkL = Uzdz{detK)^^Tr [A-D,A- + A-^^.A-^ (2.4) 

Our last step relates the reference connection and the reference bilinear structure in order to 
obtain a result which is independent upon the trivialization we used. Define the field components 
(p'' = iyiA £V , where = and plug it in (|2.4j) . One gets 

1 



%ed = -dzdzTr 



where K^j are the components of the inverse bilinear structure, that is K^jK^^ = Sl. In order to 
have a result which is independent on the trivialization, just set the reference connection to be 
the generalized Chern connection of the bilinear structure K, that is T^-. = K-jidzK^^ . Therefore, 
choosing our reference trivialization (T, K) data to satisfy this natural condition, we get 

C-red = ^dzdzTr [eijip'Dgif^] (2.6) 
which is a well defined (1, l)-form on S. Hence the action for the reduced theory is given by 

Sred = — I J^red = —\ ( dzdzTr [eijip" D^if'] . 
9s Jy, 9s ^ Jt. 

2.2 Deformations of the complex structure 

Let us now discuss certain variations of the complex structures of the manifold Cy(S, V) following 
the approach of Kodaira. 

A general complex structure variation of the vector bundle structure (|2.1jl is given by the 
deformed patching conditions 

C{a) = f{am +^(a)(/3) , W(/3)) 

= M;^am ikf^)) [^U + ik^h^iP))] (2-7) 

where a and /3 label the two local charts, and 6 and ^* are analytic functions on double patch 
intersections. The variation is trivial if it can be re-absorbed via an analytic change of coordinates. 
Notice that in the general case, the deformation functions are constrained by the chain rules of 
multiple patch intersection. 

In the following, we will consider variations leaving invariant the complex structure on S. It 
is obvious that this coincides with the general case if the moduli space of complex structures of 
S is a point. Then, from now on, we will restrict to variations of the form 

Z{a) = f{am{z{f3)) 



6 



(2.i 



Notice that the deformed complex structure preserves the CY condition if in any U(^a) ^ ^{0) 
have det (1 + d^) = 1, where (1 + d'^fj = 5] + dj^\ 

The solution of the above CY condition can be easily given in terms of a set of potential 
functions (one for each double patch intersection modulo triple intersections identities) which 
generates the deformation, as 

where we defined = ^{a) + (^^'^^ '^{/3))- 

For later application, let us specify the previous general construction for S = P^. The 
patching on the sphere allows a drastic simplification of the above formulas. In this case the 
moduli space of complex structure of the base Riemann surface is point-like, so keeping it fixed 
is not a constraint. The sphere can be described by the standard charts Us/j\i around the north 
and south poles and the single intersection Us CiUn = is the cylinder C \ {0}. 

As it is well known, by Grothendieck's theorem, any holomorphic vector bundle on can be 
presented as a direct sum of line bundles. In our case therefore V = 0{—ni) © 0{—n2), where 
we denote by 0{—n) the line bundle defined by the glueing rules 

zn = —Zg^ and wn = z^ws- 

The CY condition for the total space Cy(S, V) is therefore ni + n2 = 2. 
The generic variation is 

a;iv = 4'[o;K^M^5,u;5)] (2.9) 

where ujs = (w^, w^). In Eq. (|2.9|) . since there is one single double intersection and no triple ones, 
the functions are just constrained to be analytic on x C^, that is are allowed to have poles 
of finite orders at and oo and have to be analytic in ujs- The relevant terms in i.e. the 
ones representing true variations of the complex structure of the initial space, are the ones which 
cannot be re-absorbed by analytic reparametrizations of los and u)n- Moreover, the Calabi-Yau 
condition in the new complex structure is solved by a single potential function X = X{zs,ujs) 
such that 

eijw'^dw^ = eijiv'^duj^ - dX, (2.10) 

where, as before, w"^ = io"^ + ^'*(z5,a;). 

Let us note now that we can, as it is usually done in Kodaira-Spencer theory, relate the 
deformed and the original complex structures by a singular change of coordinates. For the case 
at hand, it is enough to do it along the fibers above the south pole patch, namely 

w}^ = uj%, and w'g = uji. + (,zs,ujs) (2.11) 

In the singular coordinates {z,w^) the patching rule is the original linear one. Therefore, Eq. 
(|2.11|) defines naturally the transformation rule for generic sections in the deformed complex 
structure from the singular to the non singular coordinate system. 

2.3 Reduction over S = in the deformed case 

Let us now perform the reduction to the brane of the open string field theory action on a Calabi- 
Yau deformation of CY{^,V) with T. = and V = 0{-ni) © 0{-n2) with m + n2 = 2. 
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Actually, from the perspective we adopted so far it turns out that performing it is not crucially 
different from the linear case. That's because we can proceed by performing the reduction of 
the hCS theory in the singular coordinates 1)2. following the prescription proper to the linear 
undeformed case and then implement the variation of the complex structure by passing to the 
non singular variables by the proper field redefinition. 

Let us start for simplicity with the reduction in the Abelian U{1) case. In this case the cubic 
term in the hCS Lagrangian is absent and the reduction is almost straightforward. In the singular 
coordinates we obtain that 

r 



1 



jLf'dsf^dzdz (2.12) 

in both the north and south charts. The coordinate change for the fields 93* in terms of the ones 
corresponding to the deformed complex structure is induced by p.llj) . Let us recall that the 
functions ^I'* defining the deformation are built out from the potential X as in eq. (|2.1U|) . This 
expresses exactly our Lagrangian terms (patch by patch) 



(2.13) 



where Xj^ = and Xg is an arbitrary analytic function of the 0's in and of z in ((^* are 
akin to the coordinate singular coordinate of the previous subsection, while (j) stem from w*). 

The above potential term X gives the deformation of the action due to the deformation of 
the complex structure. Specifically, we have that 



Sred 



1 



/ Xs{C-Ted)s + / XAf('C 
JUs JUn 



redJN 



(2.14) 



where we explicitly indicated the resolution of the unity on the sphere 1 = Xs~^Xn- For simplicity 
we choose the x's to be simply step functions on the two hemispheres. Substituting 1)2. 12(1 and 
(fTT^ in (I^TH) we then obtain 



red 



2gs 



eijcl)^dzcl>^ dzdz — / dzX{z,(f))dzdz 



D 



(2.15) 



where D is the unit disk (south hemisphere). The disk integral can be reduced by the Stokes 
theorem, leaving finally 



'red 



1 

2^ 



eij(t)^dz(t>' dzdz + ® X{z,(t))dz 



(2.16) 



where ^ is a contour integral along the equator (we understand the factor l/2pii). 

Therefore, we see that the reduced theory gives a b-c {(3 — 7) system on the two hemispheres 
with a junction interaction along the equator and the identifications ()2.9() on the fields. 

The non-Abelian case is a bit more complicated than the Abelian one because of the tensoring 
with the (trivial) gauge bundle. This promotes the vector bundle sections to matrices and 
therefore unambiguously defining the potential function X in the general case is not immediate. 
In the following we show where the difficulty arises and what further constraint to the deformation 
of the complex structure is needed in order to suitably deal with the non-Abelian case. 

To see this let us perform the reduction on of the non Abelian hCS (|2.2j) . as we did in the 
Abelian case. Let us work in the singular coordinates and obtain again the pullback Lagrangian 
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we got in the linear case. Now, in order to pass to the non singular coordinates we have to 
promote to a matrix equation the change of variables ()2.13|) . This can be done by specifying 
a prescription for matrix ordering. Suppose we choose a specific ordering and denote it by P. 
Then our change of variable is 

Tr[eij^'d,ip^] = Tr[eij^%^] - d^TrX^ (2.17) 

while the cubic term gives 



Tr 



A^e^jicP' + ^'^'W + (2.18) 



It appears immediately that our result is complicated and seems to depend quite non-trivially 
on the matrix ordering prescription. Otherwise, it is well defined. 

The easiest way not to have to do with matrix ordering prescription is to have to do with 
only one matrix. Henceforth we restrict to the case in which X{z,uj) does not depend, say, on 
U2 and we proceed further. 

In this case the deformation formulas simplify considerably. Eq. (|2.10|) is solved by ^'^ = 
and ^'^ is determined by the potential by 



This condition can be written also as 

As far as the reduction is concerned, ea. H2.17|) is unchanged, since we do not need any pre- 
scription P] while ea. (|2.18|) simplify to 

Tr [A,e,0 + ^*)(0^' + ¥)] = Tr {A, [<P\<p']) (2.20) 

(where we used [(/>^, ^'^J = 0) which, as in the linear case, is the contribution needed to complete 
the covariant derivative. The last operation to obtain our final result is an integration by part 
in the derivative term. As 

in both the north and south charts, from the last term we get an additional contribution to the 
equator contour integral, that is fTrcj)^"^^. Adding it to the previously found term we get 
j^l §Tr{X + (p^^^^). This, by (l2l^ can be written just as §TrB, where d^B = ^'2. 

Therefore, summarizing, we find that in the non Abelian case on the Riemann sphere we are 
able to treat the deformations of the type 

o;^ = z^"4, and u;^ = 4+'^[a;| + 9^iB(z5,4)] (2.21) 

which corresponds to the choice ni = —n. This geometry has been introduced (in the matrix 
planar limit) by j7|. These geometries are CY for any potential B analytic in x C. The 
relevant reduced theory action is given by 



^ -1 

'■-'red 

9s 



-TT{(tP'D^<i^)dzdz + jlh:B{z,(j^)dz 



(2.22) 
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2.4 Reduction to matrix models 



For completeness, let us generalize to our case the argument of [S] to show the reduction of the 
previous action to the matrix models. In calculating the partition function of the open strings 
attached to the D-branes^ we can easily integrate out the gauge connection which implies the 
constraint = and then cf)"^ which implies the constrain dz4>^ = 0. As a result''' we get 

that the partition function 



■'red 



L'[^2,(/)\02]g-5,,d 0, [ D[(l)^]5ids(l)^)e-^ ^ ^ (2.23) 



The condition such that the equation dz(j)^ = admits solutions is n > in ()2.21|) . In such 
a case we have n + 1 independent solutions which in the south patch are the linear span of 

{z]^}i=0...n- 

Therefore, expanding = Y17=o-^i-^s l|2.23)) . we are left with the multi-matrix model 
partition function 

Zred = J n^^^e-^^^O'-^") (2.24) 

i 

where the potential is given by 

W{Xo,...Xn) = jidzB(^z,f2Xiz'^ (2.25) 

This coincides with the one introduced by [7j in the matrix planar limit. 

The original Dijkgraaf-Vafa case is reproduced for n = 0. Then, the only non trivial complex 
structure deformation in ()2.21|) is with B = ^W{u)i) (since any other dependence in z can be re- 
absorbed by analytic reparametrizations) and hence we get the one matrix model with potential 
W. 

The above formula can be also inferred by just generalizing another CFT argument by 
Dijkgraaf-Vafa to the geometry ()2.21|) . To this end let us consider again the two dimensional 
theory defined by the action 

S = - [ Tr((/>2L',-0i) (2.26) 

where D2 = dz + This is a gauged chiral conformal field theory: a gauged b-c (P-j) 

system in which (pi and (p2 are conformal fields of dimensions — n/2 and 1 + n/2 respectively. For 
any n, the fields (pi and (/>2 are canonically conjugated and on the plane they satisfy the usual 
OPE 

Mz)Mw) ~ (2.27) 
z — w 

In Hamiltonian formalism, that is in the radial quantization of the CFT, the partition function 
is given as 

Z = (out|in) . (2.28) 



^Notice that the field redefinition ip ^ (j) has unit Jacobian because of the CY condition. 

"''See for more details about how the ghost contribution in the maximal Abelian gauge compensates the det~^ 
term. 
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The deformed transformation 



</)'2 = z"+2(</)2 + 9^,5(z, </.!)) (2.29) 
is given on the cyhnder z = e*" as 

and is implemented by the operator 

U = exp(TT^^Biz,<Pi)^ (2.31) 

Therefore the new partition function is 

Z = (out|C/|in) . (2.32) 

which is our result. 

We remark that this is an a posteriori argument, it is a consistency check but does not explain 
the dynamical origin of the matrix model from the string theory describing the brane dynamics. 



3 Engineering matrix models 

Once the link between D-brane configurations and multi-matrix models is established, the next 
natural question to ask is what kind of matrix models we get in this way. In this section we 
single out what is the most general type of multi-matrix model we can engineer by deforming 
D-branes on 2-cycles in the above way and we produce some examples. 

The geometric potential B{z,uj) is a general holomorphic function on C* x C but the terms 
actually contributing to a change in the complex structure and giving a non zero matrix potential 
are of the form 

oo d-n 

B{z,u;) = Y,Y.i'^'''''^' (3.33) 

d=l k=0 

(k) 

where are the times of the potential and to is to be identified with the coordinate loi of the 
previous section. It can be easily proven that other terms in the expansion can be re-absorbed 
by an analytic change of coordinates in the geometry. Consistently with the geometric theory of 
deformations they do not contribute to the matrix potential. 

The degree of the potential B is the maximum d such that t^^^ is non-zero for some k in ()3.33|) 
and corresponds to the degree of the matrix potential, obtained as 

WiXo, ...,Xn)= j^B i^, J2 j • (3.34) 

Since this operation is linear, from H3.33() and ()3.34p one gets a matrix potential of the form 

oo d-n 

W{Xo, . . . , X„) = ^ ^ )(Xo, ...,Xn) (3.35) 

d=0 k=0 
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where each term 



wf\Xo,...,Xn)= Yl Xn---Xi, (3.36) 



corresponds to b'^\z,uj) = z~^~^ijj'^ for < d < +00 and < /c < d • n. Note that these 
are directly obtained in completely symmetric ordered form with respect to the indices ii, . . .i^ 
labeling the different matrix variables. In the following we will sometimes write simply the 
polynomial W for c-number variables W{xo, . . . ,Xn), understanding the total symmetrization 
when matrices are plugged in. 

As it was already anticipated in the previous section, the one matrix models corresponds 
directly to the Dijkgraaf-Vafa case with n = and therefore B = ^W{ui). 

Two matrix models are obtained by considering the case n = 1. Some of them have been 
derived in [7j. In this case, it is possible to engineer a general function for two commuting 
variables. In fact 



and the matrix potential reads 

W{x,,x,) = Y,Y.'f[l)x\xt-' (3.38) 

d=l k=Q ^ ^ 

which is, upon varying the possible couplings, a generic analytic potential in the two variables 
xq and xi. The only constraint is the matrix ordering which is always the symmetric one. In 
particular, it is easy to engineer a two matrix model with bilinear coupling. This is achieved by 
choosing, for n = 1, the geometric potential to be 

B{z, uj) = - \v{uj) + U (-]] + (3.39) 
z I \ z J i Zz^ 

which generates the matrix potential 

W{xo, xi) = V{xq) + U{xi) + cxoxi (3.40) 

In general, the multi-matrix models one can engineer are not of arbitrary form. Actually, on 
top of the fact that we can generate only matrix potentials with symmetric ordering, there are 
also constraints between possibly different couplings. This can be inferred from the fact that 
a polynomial function in n + 1 variables of total maximal degree d is specified by many more 
coefficients than the ones we have at our disposal. (As an example, if n = 3 and d = 3 we would 
need 10 coefficients, while we have only 4 at our disposal.) 

To end this section, we would like to remark that some deformations can connect cases 
with different values of n. The geometric equivalence of seemingly different complex structures 
becomes in fact explicit at the matrix model level. As an example, let us consider the case n = 2 
and a geometric potential of the form 

B{z, Lo) = -\z-^^'^ + z-^F{oj). 
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Out of this, one obtains 

W{xo,xi,X2) = (F'(xo) - xi) X2 + ^ iF"{xo)xl) 

After integration of X2, which appears hnearly, this theory is equivalent to a one-matrix model 
with potential 

Vixo) = If"F'\xo), 

which is equivalent to n = and B{z, u>) = jV{oj). As a matter of fact, the geometry with n = 2 
and B = —^z~*uj'^ is equal, upon diagonalization, to the geometry n = and B = 0. 



4 General properties of two— matrix models. 

In the second part of the paper we concentrate on a subclass of matrix models: the two-matrix 
models with bilinear coupling between the two matrices but arbitrary self-coupling of each ma- 
trix. Our purpose is to find exact quantum solutions. For this reason we will solve them with 
the method of orthogonal polynomials. This method allows one to explicitly perform the path 
integration, so that one is left with quantum equations. The two basic ingredients are the quan- 
tum equations of motion and the integrable linear systems. The latter in particular uncover the 
integrable nature of two-matrix models, which stems from the Toda lattice hierarchy under- 
lying all of them. Our approach for solving two-matrix models consists in solving the quantum 
equations of motion and, then, using the recursiveness intrinsic to integrability (the flow equa- 
tions), finding explicit expressions for the correlators. An alternative method is based on the 
W constraints on the functional integral. We do not use it here, but one can find definitions, 
applications and comparisons with the other methods in j^H I52| I5,S| I54j . 

For general reviews on matrix models applied to string theory, see |2U1 1211 122| . For general 
application of matrix models, see j2,S| I24j . Early literature on two-matrix models is contained 
in refs. j25j through j38j. The method used in the present paper, although implicit in the early 
literature, is, quite incomprehensibly, seldom utilized. Different methods (from the saddle point 
to loop equations) are often preferred, see refs.jSH] through [SO]. 



4.1 Review of old formulas 

The model of two Hermitian N x N matrices Mi and M2 with bilinear coupling, see ()3.39I3.4U() . 
is introduced in terms of the partition function 

Z^it, c) = j dAfidMae*^^, W = Vi + V2 + CM1M2 (4.1) 

with potentials 

Pa 

K. = 5^WM; a = 1,2. (4.2) 



r=l 



where pa are finite numbers. These potentials define the model. We denote by M.p^^p2 the 
corresponding two-matrix model. 
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We are interested in computing correlation functions (CF's) of the operators 



Tk = trM^ Gk = trM^, V/c, 
For this reason we complete the above model by replacing H4.2() with the more general potentials 

oo 

Va=^ta,rM'a, a = 1,2 (4.3) 

r=l 

where ta,r = ia,r for r < p^. The CF's are defined by 

Qn+m 

< r., . . . . . . >= . . . . . . ^) (4-4) 

where, in the RHS, all the t^^j. are set equal to ia,r for f ^ Pa and the remaining are set to 
zero. The unusual + sign at the exponent of the integrand in 1)4. II) is because we want to use a 
uniform notation for physical couplings ta^r and sources ta^r (for the convergence of the integrals, 
see below). From now on we will not distinguish between ta,r and ta^r and use throughout only 
ta,r- We hope the context will always make clear what we are referring to. 

We recall that the ordinary procedure to calculate the partition function consists of three 
steps , [SH] , |SZ| : {i) one integrates out the angular part so that only the integrations over the 
eigenvalues are left; (ii) one introduces the orthogonal monic polynomials 

Cn(Ai) = A" + lower powers, ^?n(A2) = A2 + lower powers 

which satisfy the orthogonality relations 

dAidA2en(Ai)e^i(^i)+^2(^2)+^^i^27y™(A2) = hn{t,c)6nm (4.5) 



{iii), using the orthogonality relation 1)4.51) and the properties of the Vandermonde determinants, 
one can easily calculate the partition function 



N-l 



Z^it, c) = const iV! Jl hi (4.6) 

i=0 

whereby we see that knowing the partition function means knowing the coefficients hn{t, c). 

The crucial point is that the information concerning the latter can be encoded in 1) a suitable 
linear system subject to certain 2) equations of motion (coupling conditions), together with 3) 
relations that allows us to reconstruct Z^r. 

Let us introduce some convenient notations. We will meet infinite matrices Mjj with < 
i,j < 00. For any such matrix M, we define 

M = H-^MH, Hij = hi6ij, Mij = Mji, Mi{j) = Mjj_i. 

We represent such matrices in the lower right quadrant of the (i, j) plane. They all have a band 
structure, with nonzero elements belonging to a band of lines parallel to the main descending 
diagonal. We will write M G [a,b], if all its non-zero lines are between the a-th and the 6— th 
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ones, setting a = for the main diagonal. Moreover M_|_ will denote the upper triangular part 
of M (including the main diagonal) , while -/Vf_ = M — M4. . We will write 



Tr(M) = ^ M, 



T\A.. 
II 



j=0 

Next we pass from the basis of orthogonal polynomials to the basis of orthogonal functions 
*n(Ai) = e^^(^^)C„(Ai), ^ni\2) = e^^(^^)r?n(A2). 
The orthogonality relation ()4.5() becomes 

y d\rd\2^n(\x)e'^'^^^mi\2) = SnnMt, c). (4.7) 

We will denote by ^ the semi-infinite column vector (^'q, ^i, ^'2, • • • , )* and by $ the vector 
($0) ^1, ^2, •••,)*• Then we introduce the following Q-type matrices 

j dXidX2^niXl)Xae''^'^^<i>miX2) = Qnmia)h^ = Qmnia)hn, a = 1,2. (4.8) 

Beside these Q matrices, we will need two P-type matrices, defined by 

I dXidX2[^^-i'n{Xl))e''^'^'<^m{X2) = Pnm{l)hm (4.9) 

dAidA2%,(Ai)e=^i^2 (^_^,„(A2)) = Pmn{2)K (4.10) 



For later use we also introduce 

d 



dAidA2(^en(Ai))e^i(^^)+^^(^^)+^^i^^r/^(A2) ^ /'„°™(l)/im (4.11) 



dAi(iA2e„(Ai)e^^(^^)+^^(^^)+^^^^^(^7?^(A2))=P^j2)/i„ (4.12) 
Let us come now to the three elements announced above. 

1) Quantum equations of motion. The two matrices 1)4. 8p we introduced above are not in- 
dependent. More precisely both Q(a)'s can be expressed in terms of only one of them and one 
matrix P. Expressing the trivial fact that the integral of the total derivative of the integrand 
in eq. (|4.7j) with respect to Ai and A2 vanishes, we can easily derive the constraints or coupling 
conditions, or quantum equations of motion 

P°(l) + Vi + cQ{2) = 0, cQ(l) + Vi + V°{2) = 0, (4.13) 

These may be considered the quantum analog of the classical equations of motion. The difference 
with the classical equations of motion of the original matrix model is that, instead of the N x N 
matrices Mi and M2, here we have infinite Q{1) and Q{2) matrices together with the quantum 
deformation terms given by P°{1) and V°{2), respectively. From the coupling conditions it follows 
at once that 

Q{a) e [-ma,na], a = 1,2 
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where 



mi = p2 



1, 



m2 



n2 = pi 



1 



where pa, a = ^, 2 is the highest order of the interacting part of the potential Va (see (|4.2() ). 

2) The associated linear systems. The derivation of the hnear systems associated to our 
matrix model is very simple. We take the derivatives of eas. H4.7|) with respect to the time 
parameters t^^r, and use eas. (|4.8() . We get in this way the time evolution of ^ and which can 
be represented in two different ways: 



Discrete Linear System I: 



( g(l)^'(Ai) = Ai*(Ai), 
^^(Ai) = Qni)+^(Ai), 



d 



■^(Ai) = -Q'=(2)_vI/(Ai), 



V |^^i/(Ai) = p(i)^(Ai: 



The corresponding consistency conditions are 



^ Q(1) = [Q(1), QHa)-] 



a 



a,k 



1,2 



(4.14) 



(4.15a) 
(4.15b) 



In a similar way we can get the time evolution of ^ via a discrete linear system II, whose 
consistency conditions are: 



[Q(2), P(2)] = l, (4.16a) 
Q{2) = [Q'ia)+, Q(2)] (4.16b) 



dta,k 

One can write down flows for P{1) and P{2) as well, but we will not need them in the sequel. 

3) Reconstruction f ormulae. The third element announced above is the link between the 
quantities that appear in the linear system and in the quantum equations of motion with the 
original partition function. We have 

lnZNit,c) = Tr(Q''{a)), a = 1,2 (4.17) 



dta 

It is evident that, by using the equations H4.15bl4.16b|) above we can express all the derivatives 
of Zn in terms of the elements of the Q matrices. For example 

-lnZ7v(t,c)= Q'X") , a = 1,2 (4.18) 



dti^ldta,r ^ J N,N-l 

and so on. We recall that the derivatives of F{N,t,c) = lnZAr(t, c) at prescribed values of the 
coupling are nothing but the correlation functions of the model. 

The above derivation is rigorous when, for example, highest negative even couplings guarantee 
that the measure in ()4.5() is square-integrable and decreases more then polynomially at infinity. 
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But for generic values of the couplings it is heuristic. Nevertheless we notice that the consistency 
and quantum equations of motion make sense for any value of the couplings, and also when the 
couplings are infinite in number. In the latter case eas. H4.15b|) and eas. (|4.16b|) form nothing but 
a very well-known discrete integrable hierarchy, the Toda lattice hierarchy, see j58j . 

From these considerations it is clearly very convenient to refer to the integrable system for- 
mulation rather then to the original path integral formulation of our problem. This allows us not 
only to extend our problem to a larger region of the parameter space, but also to make full use 
of integrability. 

To end this section, we collect a few formulas we will need later on. First, we will be using 
the following coordinatization of the Jacobi matrices 



mi m2 



Q(l)=/+ + ^^a;(i)^M-^ Q{2) = 1+ + Y,Y.^i{i)Ei,i-i (4.19) 

i 1=0 i 1=0 

where 1+ = X^i=o iEij)k,i = 5i,k^j,i- One can immediately see that 

(g+(l)) + ao(i)<^ij, {q_{2))^= R{i)5j,i.i (4.20) 

where R{i + 1) = As a consequence of this coordinatization, eq. (|4.18|l gives in particular 

the two important relations 

F{N,t,c) = ai{N), (4.21) 



and 

92 



dti.idt2: 



-F{N,t,c) = R{N) (4.22) 



For reasons of brevity we do not even touch on the subject of W-constraints. The latter 
are constraints on the partition function under the form of algebraic structure, see |51[ I54j for 
instance. They are obtained by putting together quantum equations of motion and flow equations. 
W-constraints (which are also called loop equations or Schwinger-Dyson equations) can be used 
to solve matrix models, but such a procedure is less efficient than the one used in the sequel. 



4.2 Homogeneity and genus expansion. 

The CF's we compute are genus expanded. The genus expansion is strictly connected with the 
homogeneity properties of the CF's. The contribution pertinent to any genus is a homogeneous 
function of the couplings (and N) with respect to appropriate degrees assigned to all the involved 
quantities. Precisely, we assign to the couplings the following degrees 

deg{) = [l = x(l - fc), [N]=x, [c] = -x (4.23) 

where x is an arbitrary positive number. Here and in the following is treated as a coupling 
^1,0 = t2,o- If we rescale the couplings as follows 
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we expect the free energy to scale like 

oo 

F^^A^'(2-2/')F/, (4.24) 

h=Q 

where h is the genus. In other words 

[Fh] = 2x{l - h) (4.25) 

Fh is interpreted as the result of summing the open string partition function at fixed h over all 
the boundaries. 

The CF's will be expanded accordingly. Such expectation, based on a path integral analysis, 
remains true in our setup due to the fact that the homogeneity properties carry over to the 
Toda lattice hierarchy. To this end we have simply to consider a genus expansion for all the 
coordinate fields that appear in (5(1) and (5(2), see (|4.19[ I4.2()|l . The Toda lattice hierarchy 
splits accordingly. In genus the following assignments 

[a}°)] = (/ + l)x, [hf^] = {l + l)x, = 2 X (4.26) 

correspond exactly to the assignments 1)4. 23(1 and [Fq] = 2x. 

5 Solving two— matrix models 

As already pointed out in the previous section, a way to solve a two-matrix model is to solve 
the coupling conditions (quantum equations of motion). This allows us to determine the 'fields' 
ai{n), bi{n) and R{n). Once these are known we can compute all the correlation functions starting 
from (j4.17|) by repeated use of eqs. (|4.15al4.15bl4.16a|4.16b|) . which form the flows of the Toda 
lattice hierarchy. As for the free energy F(t,N,c), see for instance f j52jl. In using this 
method, the bi-Gaussian model A^2,2 was solved. This is of course a simple model. However it 
is useful to check the coincidence of the results obtained in this way in the decoupling c = case 
with the available results obtained by the traditional method based on eigenvalue density and 
resolvent for the Gaussian one-matrix model. 

A very interesting case is the model Mofi, i.e. the limiting model when only the c parameter 
is different from zero. As a path integral this model does not make much sense. However, as we 
saw above, it does have sense as integrable system to which the appropriate coupling conditions 
are applied. It turns out that this model describes c = 1 string theory at the self-dual radius, as 
was shown in |521 153j . 

5.1 Solving the quantum EoM's: A4^^2 model 

The next model in order of complexity is the A^3,2 model, [SH]. The relevant quantum equations 
of motion are 

P°(l) + 3hQ{lf + 2t2Q(l) +ti + cQ{2) = (5.1) 
P°(2) + 2s2(5(2) + si + cQ(l) = (5.2) 
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Using the coordinatization (|4.19|) and H4.20() they produce the followmg equations for the fields 
ai{n),bi{n),R{n) 



c62(n) + 3t3R{n)R{n - 1) = 

2t2R{n) + c6i(n) + ?,UR{n) {ao{n) + ao(n - 1)^ = 

3i3 {aQ{nf + ai(n) + ai(n + 1)) + 2t2ao(n) + ti+ cbo{n) = 

n + 3t3ai(n) ^ao(?T.) + ao{n — 1)^ + 2t2ai(n) + cR{n) = 
2s2R{n) + cai(n) = 
2s2&o('^) + si + cao(n) = 
n + 2s2fei(n) +cR{n) = 



(5.3) 



where we have introduced the simphfied notation 



t2,k = Sk 



One easily realizes that the second, fourth, fifth and seventh equations are linearly dependent. 
The remaining equations determine the lattice fields uq, ai,bQ,bi,b2, R completely. The M.3^2 
model, even though comparatively simple is already rather complex due to the large number of 
involved fields. Therefore, for pedagogical purposes, let us further simplify it, by setting c = 
and considering only the ao,ai fields. This corresponds to the one-matrix model with cubic 
interaction. In eas. ()5.3(l . for consistency, we have to set also R = 0. The relevant equations are 



One can derive ai from the second equation in terms of ao and replace it into the first. One 
gets in this way a cubic algebraic recursive equation for oq. We solve it with a genus by genus 
approach. The first step is to start with genus 0. To do so one simply ignores the increments ±1 
on the n entry. In this way we get an ordinary cubic algebraic equation in the unknown ao: 



where, for simplicity and without loss of generality, we have set ti = 0. In the large N limit, it is 
convenient to shift to the continuous formalism, by defining the continuous variable x = ^. It is 
also useful to make contact with section 4 of for a comparison. So, also in order to simplify a 
bit further the notation, we set ^2 = — y and ^3 = —Ng, where g is the cubic coupling constant 
there. Moreover we denote / = 3gao. Then ea. ()5.5() becomes 




(5.4) 




(5.5) 



(5.6) 



18g^x + f{l + f){l + 2f) = 



(5.7) 



It is easy to find the three solutions, which for small x take the form 



/i = -18g'^x - 972g^x^ - 933125^x3 - 11022480/x^ + 0{x^) 
/2 = -1 - I8g^x + 972/x2 - 93312/x^ + 11022480/x^ + O(x^) 

/3 = + mg^x + 186624/x^ + O(x^) 



(5.8) 
(5.9) 

(5.10) 
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From them we can easily write the fields qq and ai in terms of g and x. They therefore lead to 
three different solutions for the correlators. Later on we will show how to compute the latter. 
But, before, let us discuss the meaning of these three solutions. To start with, comparing this 
with inS]) we see that the first solution corresponds to the unique solution found there provided we 
set X = 1. It corresponds to the minimum of the classical potential. In fact the correspondence 
with fS^* can be made very precise: one can easily verify that eqs.(46) there are nothing but 
eas. (|5.5l5.6() provided we make the identifications: a + b = and (6 — a)^ = 4g(Oi with x = 1. 
In the interval (2a, 26) represents the cut in the eigenvalue distribution function. This cut 
therefore has to be found in our formalism in the (ao,ai) plane, once we forget the dependence 
of the latter on x. 

The classical potential for the continuous eigenvalue function A(x) (which is A,„/\/iV in the 
large limit), is Vd = + gX^ . It has extrema at A = and X = —l/3g. To find the classical 
limit we have to drop the last term in ea. H5.5|) . Remember that dropping the latter (with c = 0) 
is equivalent to writing ea. H5.1|l as V({Q{1)) = 0. In other words the latter is the equation 
that identifies the extrema of the quantum potential Vi. They are three, / = 0,-1,— 1/2, 
which corresponds to oq = 0, —l/3g, —l/6g, not two as in the classical case. So we see that bq 
approaches, in the classical limit, the classical eigenvalue function. Moreover / = corresponds 
to the minimum of the potential, / = — 1 corresponds to the maximum and / = —1/2 to the 
vanishing of the second derivative. The first two cases (the classical extrema) are characterized 
by the fact that oi = 0, while the third corresponds to a non- vanishing oi. 

From this simple example we learn three important pieces of informations. 

• The number of solutions of the quantum problem (i.e. the number of solutions to ea. (|5.7() ) 
is larger than the number of the extrema of the classical potential. 

• The field oq can be regarded as the quantum version of the classical eigenvalue function. 

• The classical extrema are obtained by setting, together with n = in ea. H5.5|l . also ai = 0. 

As we shall see, the last condition, in the most general case, must be replaced by all fields 
ai,a2, ... present in the problem being set to zero (except oq). 

What is the meaning of the third solution, / = —1/2? It is a non-perturbative solution. It 
cannot be seen in the saddle point approximation. More on this later on. 

It is now easy to extend the analysis to the full A^3,2 model. In genus (i.e. disregarding the 
±1 increments on n) it leads to the following set of equations, see j54j : 

ai(n) = 
bi{n) = 

and the recursion relations 

2ao(n) 
2R{n) 

As expected the last two equations lead to the same cubic equation for oq as (|5.7|) with modified 
coefficients. It is interesting to find the classical extrema of this model. According to the above 



-^i?(n), 
c 

n + cR{n) 

2^2 



6o(n) = 
62 (n) 



coo(n) 
3i3 



R{n) 



(5.11) 



2t2 ^ c 
3i3 6^2*3 



c + 



n 



R{n) 



anin) + 7^ — 

2S2 V6S2t3 125^*3 



ao(n) 



(5.12) 
(5.13) 
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recipe we must set ai = 61 = 62 = 0, and drop the first term in the third and sixth of eas. H5.3|) . 
In the genus version of the latter this leads to i? = and to 

3*3 ag + 2*2 ao + cbo = 

2s2bo + cao = (5.14) 

The extrema correspond therefore to 

— 4^2*2 

a = 0, and ao = (5.15) 

652*3 

which is what one expects by completing the quadratures in the original classical potential. Of 
course, as above, in this way we find only two extrema out of three. 

The solutions of the 7^3,2 model therefore can be found in the same way as in the simplified 
one-matrix model above. 

5.2 Solving the quantum EoM's: A^4 2 model 

As a further example we briefly analyse the ^^4,2 model with, for simplicity, only the t4, t2,c and 
S2 couplings switched on. The quantum equations of motion are 

P°(l) + 4UQ{lf + 2i2Q(l) + cQ{2) = (5.16) 
P°(2) + 2s2Q(2) + cQ(l) = (5.17) 

The matrices Q{1) € [-1,1] and Q{2) G [-1,3]. Using the coordinatization and 

they produce the following equations for the fields a/(n), 6/(n), R{n), which we write down in the 

genus version: 

n + 12*4(0? + agai) + 2*201 + cR = 
4*4(00 + 600O1) + 2*200 + c6o = 

12*4(0^ + 01) + 2*2 + = 

12*400 + C— r = 

n + 252^*1 +cR = 
2S2&0 + coo = 

2s2i? + cai = (5.18) 

Now let us proceed as in the 7^3,2 and set c = 0. We obtain two decoupled one-matrix model 
systems, a Gaussian one on the right and a quartic one on the left. We are interested in the 
latter. The relevant equations are 

n + 12*4(01 + OoOi) + 2*201 = 

4*4(0^ + 60001) + 2*200 = (5.19) 
Now, the second equation admit the solution oq = 0. Replacing it into the first we obtain 

n+ 12*401 + +2*201 = (5.20) 
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If 7^ 0, we can derive Oi from the second equation and replace the result into the first, obtaining 
a biquadratic equation for ao 

12^^-20(|i)'a^-12^ag-l = (5.21) 

Both (|5.2Uj) and (|5.21j) can be solved exactly. They give rise to six distinct (in general complex) 
solutions. The classical potential for the eigenvalue function in the large N limit is V4(A) ~ 
t2A^ + t^X^. This potential has one or three real solutions depending on whether t2 and t4 have 
the same or opposite sign: A = and ^^'^ = In order to single out among the above six the 

solutions that correspond to those in the classical limit, we follow the above given recipe. We 
drop the first term in the first ea. ()5.19|) and set ai = in both. We are left with 

Atiaff + 2t2ao = (5.22) 

This gives exactly the expected classical extrema for oq. Once these are determined we can easily 
find the corresponding quantum solutions either in exact form or in series of x = n/N. Following 
the example of the previous subsection we can also determine the solutions of the complete 7^4,2 
model. 



5.3 Solving the quantum EoM's: A^3 3 model 

We study the model in the case ti = si = and limit ourselves to writing down the genus 
quantum equations of motion: 

St^cttQ + 2t2caQ — SGsst^boR + c^bo — 1252^3-^ = 
3s3c6o + 2s2c6o — SGsstsOo-R — 12s3t2^ + oqc^ = 

nc + Rc^ - ISsst^R^ - SGs^tsaoboR - 12s2t3aoR - 12^25360^ - 4^2*2-^ = (5.23) 

_3£3^2 

9 

3t3^2 

9 

For simplicity we compute only the classical vacua. To this end we drop the first term in the 
third equation and solve the resulting system. Then we set 02 = 62 = ai = 61 = as well as 
R = 0. In this branch we therefore have 

Stgag + 2*200 + c6o = (5.24) 
3s3bl + 2S260 + cao = (5.25) 

From the first we can get 60 = — ^(Stsag + 2*200)) whence we get either oq = or the cubic 
equation 

27s3tlal + 36*253*300 + (1253*2 - 6052*3)00 + c(c^ - 4S2*2) = 

Therefore in general we have four classical extrema, with non-vanishing oo and 60 while all the 
other fields vanish. A series expansion about these solutions is easy to find. For instance around 



6S3, „ 2S2 

oi = boR R, 02 

9 9 

bi = — ^-^oo-R -R, 62 = 

9 9 
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the vacuum uq = bo = R = we have 



(c2 - 4s2i2)2 

648[-8csl4 + c^slt2{c^ + 852*2)^3 - 4cs|s3(c2 + 252*2)^3 + ^Qs^t^jx'^ 

^ (c2 - 4s2i2)^ 

12(-2g3ti + CS2t3)x 

° (c2 - 4s2t2)2 

648[16s|f| - 4cs§t|(c2 + 2s2i2)*3 + c'^S2S3{c^ + 852*2)^3 - 8cs^tl]x^ 



3\ 



(c2 - 452^2)5 

ca: , 18c[-16s2t3 + cg3(c2^i2s2t2)^3-16git2]x2 ^ 3 

c2-4.2t2 (C2- 4.2*2)'^ ^""^ 



For reasons of space we have Hmited the expansion in x = n/N to the quadratic order. From 
(|5.23|) one can easily compute the expansions for ai, 61, 02, 62- 



5.4 The Mp^^p^ model 

In the general case the matrix rank for Q{1) and Q{2) was given in the previous section and 
the quantum EoMs become of course very complicated. It is however simple to write down the 
equations that identify the extrema with classical analog. They are 

Vliao) + cbo = 0, V^{bo) + cao = (5.26) 

while all the other fields are set to zero. We have in general (pi — l)(p2 — 1) solutions of this 
type in perfect correspondence with the classical analysis. The simplest solution is oq = 60 = 0. 
Other solutions may be hard or even impossible to determine explicitly. Anyhow, once one such 
solution is known it is possible to find explicit expressions for the fields around it in terms of 
X = n/N. 



5.5 Calculating the correlators 

Once we know the fields ai,bj,R in a given model there exists an algorithmic procedure to 
determine the correlators. This in turn is due to the integrability underlying the Toda lattice 
hierarchy. In this subsection we give a few examples of exact correlators. The general scheme 
is known, it was already presented in |5H I54j. A few explicit examples were worked out for the 
A4o,o model in [521 153j . In these references one can find explicit calculations of correlators for 
finite A'^ and for any genus. In this subsection we limit ourselves to large genus correlators. 
To start with let us briefly review the continuous versions of the quantum equations of motion 
and the flow equations in this case (which are known as the dispersionless Toda lattice hierarchy 
flows) . 

We proceed as in |51| . First we define the continuum quantities in the following way 



^ j.ren _ ren _ £fc ren _ 

AT' ~ AT' "'^ ~iV' ~ AT 



in the large A^ limit. They are the renormalized coupling constants. In the following however we 
will understand the superscript ren. Further we define 

Fo(x)= hm C= lim/. 
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where Fq is the genus zero free energy. The second equation is merely symbohc and simply means 
that C is the continuum counterpart of /+. If we define a matrix p = nEn,n, it is easy to see 
that we have 

[I+,p]=I+, (5.27) 
The continuous counterpart gives the following basic Poisson bracket 

{C,x} = C- (5.28) 
This allows us to establish the following correspondence 

N[ , ]^{ , } (5.29) 

and, similarly, 

-Tr = - ^ / dx (5.30) 

-'^ 

together with the replacements 

Q{l)^L, Q{2)^L; 

where 



L = c+E«'r^ ^ = f + E4c^- (5-31) 

Here ai and bi are the continuum fields that replace the lattice fields ai{i) and bi{i) of eq. (|4.19|) . 
We stress that the above replacements holds only in genus 0. 

Now let us come to the fiow equations: The dispersionless limit of the extended 2-dimensional 
Toda lattice integrable hierarchy ^.15b . 4-16b ) is 

j- = {L, {L^)-}, — = {L,(L'=)_}, (5.32a) 

Otk OSk 

Here the subscript + denotes the part containing non-negative powers of C) while - indicates 
the complementary part. 

Next, the continuum version of (|4.17j) provides the link between the free energy and Lax 
operators, i.e. 

J-F = £ {L%){y)dy = £ {L%){y)dy (5.33) 

where subscript '(g)' means that we select the coefficient of the zero-th power term of This 
formula opens the way to calculate CF's by simply differentiating both sides with respect to the 
appropriate coupling constants and use eas. 1)5.31^ 15. 32b1) . Therefore this equation together with 
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the integrable hierarchy and the quantum equations of motion completely determines the genus 
zero correlators. For instance 



< -Tkn >-- 



dtkdti 



< TkCTl >-- 



dtkdsi 



and 



d^F 
dskdsi 



{L')>idL^ 



(5.34) 



(5.35) 



where § represents the residue at the simple pole in C,. 

It is now very easy to extract explicit expressions for correlators in various model. Here we 
limit ourselves for simplicity to a simple example, the Mz,2 model in the decoupling limit studied 
in subsection 2.1. In this case we have two fields ao and ai and explicit expressions for them. 
One starts from eqs. (|5.8|5.9|5.1()|) . then computes oq = — ^ and finally ai from eq. (|5.Hj) . The 
Lax operator is given in this case by 



L = C + ao + 



ai 



(5.36) 



Inserting it in the previous formulas we get 

k 



< Tk >-- 



21=0 



{k-2l)\l\l\ Jo 



and 



< TkTr > = 

r;-li 

r y 2 ' k—1 n 
1=0 p=0 n=0 q=0 



(5.37) 



(5.38) 



k\ 



{k — 1 — n)!(n — g)!g 



r—l+k-n-l l—p+n-q 



■l-n,-l 



0'l52p+2q-l-n,l] 



where [a] denotes the integral part of a. Replacing eas. (|5.8l5.9l5.1Ujl . as indicated above, we find 
explicit expressions for the correlators in terms of x and g. 

Similarly one can compute the three point functions (referred to generically as Yukawa cou- 
plings) 



n EE 

i=l \k,=Op,=0 



hi 



(5.39) 







af ^^^^ ( a'r'^'af ) (k, - 2pi)(fc2 - 2p2) e{k2 - 2p2) e{2p^ - k^) 



-a, 



h+h-ki-ks P1+P3 







1 ) (^1 - 2pi)(A:3 - 2p3) eiki - 2pi) ei2p3 - k^) 



J_j2+h-k2-k3 „P2+P3 
"r"0 



af (a[,i-'^af ) {k2 - 2p2){h - 2p^) e{2p2 - fes) e{2p3 - h) 



where 9{x) means for x < and 1 otherwise. 
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5.6 Higher genus 

We would like to introduce in this subsection a few basic notions concerning higher genus corre- 
lators. They are introduced here in order to render the subsequent discussion as self-contained 
as possible. For a more complete treatment see [S^ 15^ 151] . 

Let us consider the first example above, the decoupled model, whose general solutions 

are characterized by ea. ()5.4() and genus ones are explicitly given in eas. H5.8l5.9l5.1U|) . Let us 
start from a given genus solution, specified by oq = vq and oi = sq, where r^/Sg is anyone 
of the three solutions H5.8I5.9I5.1U|1 . Then we expand the full solution in series of e = as 
follows: 



ao = 2^ e"r„(x), = 2^ e"sn(a;) (5.40) 

n=0 n=0 



Moreover /(nit 1) is replaced by e^^^"^ f{x) for any lattice function f{n). Inserting this into (|5.4|1 
we obtain the genus equations and an infinite series of relations for the higher order terms, 
which can be recursively solved. For instance, the next to leading equations are 

etsron + St^isQ + 2si) + 2t2ri = 

Qhisiro + son) - Shsor'o + 2t2Si = (5.41) 
where a prime denotes derivative with respect to x. They can be easily solved and lead to 

_ _3 3t3gor^ + ^2^0 + Stsrps'Q 
2 ^t2 + 6*2*3^0 + 9ti(r2- so) 
^ ^ 3^ ^ 3h{s'Q + ror'o)+t2r'o 
2 ^^^2 ^ 6rot2t3 + 9(r2 - so)ti 

Similarly, for the second order we get 



3 t2{2rf + s'o' + 2s[) + 3t3{2r[so - r'^^so + 2r^,si - Ansi + ro{2rj + + 25^)) 



''^ 4^ tl + 6rot2t3 + 9tl{r^-so) 

^ 3 -3sot3{2rj + sl + 2s[) + (^2 + 3rot3){r'^so - 2r[so - 2r'^si + 4riSi) 
4 ' 9sotl - {t2 + 3rot3)2 

and so on. Replacing these expressions into the appropriate formulas for the correlators, we can 
write down their explicit genus expansions. To this end one should recall that the appropriate 
expansion for correlators is given by a power series in e, as one can infer from the free energy 
expansion 

00 

F{x,e)=Y,Fn{x)e^'' (5.42) 

/i=0 

where h is the genus. 

We give, as an example, the genus 1 contribution to < >, which is the coefficient of in 
the e expansion: 



fX 

< >i = I dy 



4 J \^{-2) • ^ -^(-1) ' 0^0 ^ ^(0) 
(1) *0' O'^O ^(2) V*0'0 -^^OiQ ^^(3) *0'0*0T^^(4) *0*0 
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+ (3) {^Ufo + ^(oT'(^o' + ^0 + Hn) 

+^"-^y^(2roSo - soro - SoVq + SrQSi + 3sQri) + ^"27^('5o - sqSq - Sq + 3sqSi] 

+(2) + + + - 

+ (1) (S^^ + ^a^^ 



(5.43) 



where tq, so,ri, si,r2, S2 are function of y and a prime denotes differentiation with respect to y. 
Moreover 

^ffc) = ^^E ^ („_2p-A;)!(p + A:)!p!''o''''~''o ^^'^^^ 

6 Two— matrix models and multiple brane configura- 
tions 

All the examples of the previous sections represent, according to the geometric description of 
section 2, the physics of D-branes wrapped around the two-dimensional sphere located in 
one of the vacua. A related problem is to describe a more complex situation with Ni D-branes 
at one vacuum and N2 = N — Ni at another. There may of course be even more complicated 
configurations with several groups of D-branes in different vacua. Let us call them multiple brane 
configurations. However the example with two groups of D-branes will be sufficient to illustrate 
the salient features of the problem. Let us consider once again the 7^3,2 model in the decoupling 
limit so that we can work with explicit formulas. We refer in particular to ea. H5.7|) . which we 
rewrite here in the form 

18g'^x + z{l + z)il + 2z) = (6.45) 
This can be solved exactly for z and gives the three solutions 

1 1 I(x) 

2 4/(x) 12 ^ ^ 
1 l-iVs l + iV3I(x) , , , 



1 /3 

l{x) = 3^/^ f-324c/2x + ^3^-1 + 3499254x2') (6.49) 



where 



When expanded for small x they give rise to the series (|5.8|5.9l5.10j) . respectively. However this 
is not a very illuminating way of studying ea. 1)6. 45(1 . The best way is to consider it a plane curve, 
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j59j . in the complex z, x plane. Then it represents a genus Riemann surface with punctures at 
X = and x = oo. It is made of three sheets joined through cuts running fromz = -l/iVsiOSg"^) 
to z = l/(108\/3g^). The solutions (|5.8I5.9I5.1U|) correspond to the values z takes near x = 0, 
away from the cuts. Therefore we can pass from one solution to another by crossing the cuts. We 
call the Riemann surface so constructed the quantum Riemann surface associated to the model. 
This Riemann surface picture is the clue to understanding the solutions with multiple brane 
configurations. Let us see how. 

Let us start from A^3,2) set c = and concentrate on the cubic interaction part. The relevant 
equations are (setting ti = 0) 



3^3(ao(?^)^ + ai{n) + ai(n + l)j + 2^2^^o(?^) = 

n + 3t3ai{n)(ao{n) + ao{n - 1)) + 2t2ai(n) = (6.50) 



This equation in genus has three solutions. Let us denote by oq = tq one such genus solution 
and ai = sq the corresponding oi. Similarly, we pick another solution and we denote it ao = uq 
and ai = vq. They represent the lowest order contribution of expansions like those considered in 
section 2.6. 



i?(n)^i?(:E) = ^eV„(:E) 

n=0 
oo 

S{n) S{x) = J2 e"«n(2;) (6.51) 

n=0 

oo 

[/(n)^[/(:E) = J^e"n„(x) 

n=0 
oo 

V{n) V{x) = J2 e"'^n(^) (6-52) 

n=0 

where we have indicated the large A'" expansion. R{n), S{n) and U{n), V{n) form, separately, two 
couples of solutions of (|6.50|) . We recall that in the analogous problem formulated in the familiar 
saddle point approach one sets A'^i eigenvalues Aj in one vacuum and in another. Considering 
the analogy between the field uq in genus and the classical eigenvalues, we are led to pose 
the following problem: does there exist a solution of 1)6. 5U() that corresponds to R{n),S{n) for 
< n < A^i - 1 and to U{n),V{n) ior Ni < n < N - 1? This means that (|U3n|) must hold 
with ao{n),ai{n) replaced by R{n),S{n) ior < n < Ni — 1 and by U{n),V{n) for n > Ni, 
respectively. However in addition we have the boundary equations 



3t3[R{Ni - ly + S{Ni - 1) + V{Ni)j + 2t2R{Ni - 1) = 

A^i + 3t3ViNi)(u{Ni) + R{Ni - 1)) + 2t2F(iVi) = (6.53) 



While all the other equations are the same as in the previous section, these two equations represent 
the real novelty: they mix two different solutions. They are obviously satisfied if it makes sense 
to identify 

V{Ni) = S(iVi), R{Ni - 1) = U{Ni - 1) (6.54) 
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In the discrete formalism this is not easy to check. Therefore we shift to the continuous formahsm. 
Recahing what we did in subsection 2.5, 2.6, in the large N limit we set N2 = aNi , N = (l+a)Ni, 
Ni/N = 1/(1 + a) = p. We also define n/N = x for < n < iVi - 1 and n/N = Ni/N + (n - 
Ni)/N = P + y. In this formalism eas. (|6.54|) read 

V{(3) = lim S{/3 + e), R{(3) = lim U{/3 - e) (6.55) 

This is nothing but the statement that at x = /3 we are crossing the cut that separates the two 
solutions. In hindsight this is quite obvious: the only way to satisfy ea. H6.5U|) with two different 
solutions is to cross the cut that join the corresponding sheets in the Riemann surface introduced 
above. 

With this recognition in mind one can now set out to calculate correlators in a theory with 
two sets of A'^i and ^^^^2 D-branes in two different vacua. Leaving a more complete treatment 
for another occasion we can easily exhibit as an example the two point correlators analogous to 
(|5.38() . The correlator is given by the sum of two terms, each one equal to the right hand side 
of (|5.38|) : in the first ao,ai are replaced by ro,so evaluated in x, while in the second they are 
replaced by uq, vq evaluated in (3 + y. 

The conclusion is therefore that quantum solutions to (|5.4|) that correspond to two groups 
of D-branes do exist. Moreover, due to the structure of the quantum EoMs which we have 
explored in the previous section, it is easy to generalize this conclusion. Every model will be 
characterized by a quantum Riemann surface with cuts that separate different solutions. It is 
therefore possible to construct multiple brane configurations, as quantum solutions of the matrix 
model, by means of solutions of the quantum EoMs on different sheets that join to one another 
across the appropriate cuts. 

7 Conclusions 

In this paper we considered B-model D-branes on 2-cycles of local Calabi-Yau geometries. The 
theory describing these objects is given by the reduction to the D-brane world-volume of the open 
string theory with Dirichlet boundary conditions on it. We have described a precise dimensional 
reduction scheme for the holomorphic Chern-Simons theory, that is the B-model open string field 
theory, to the 2-cycle. This has been done for generic local CY geometries modeled around an 
arbitrary Riemann surface. In the case of the conifold geometry, i.e. when the 2-cycle is a P^, we 
have considered the relevant effective theory and found that it is given by a multi-matrix model. 
The number of matrices involved depends on the reference complex structure about which we 
calculate the coupling to the Calabi-Yau complex moduli. The multi-matrix potential is fixed 
by the complex moduli in a well defined and simple way. The various allowed couplings turn out 
to be in correspondence with the projective parameters of complex structure deformation. The 
matrix models we have obtained are of generic type if they involve one or two matrices, while we 
found relevant constraints within their parameters for more than two matrices. 

We have studied the geometric engineering of the multi-matrix models and provided both a 
general reduction scheme and some examples. Actually, some multi-matrix models are reducible 
to models involving less matrices. This phenomenon has a clear counterpart on the geometri- 
cal side, corresponding to the fact that different reference conifold complex structures can be 
connected via specific deformations. 

In the second part of the paper we have focused on two-matrix models with bilinear couplings. 
We have illustrated a general method to exactly solve these models. It consists in solving the 
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quantum equations of motion and making subsequently use of the integrable flow equations. We 
have exhibited several examples of solutions in genus as well as in higher order approximations. 
We have also discussed the relation of our method to the more popular saddle point approxima- 
tion. One of the relevant differences is that our method leads in general to more solutions than 
the saddle point one. 

Finally we have discussed the brane interpretation of our results. They represent the deforma- 
tions of the complex structure generated by the strings attached to N D-branes wrapped around 
a 2-cycle in a local Calabi-Yau geometry. However it is possible to obtain more general solutions: 
group of D-branes localized near different vacua. This is due to a remarkable property of our 
approach: all the data of all the solutions of a given model are encoded in a Riemann surface (a 
plane curve), which we call the quantum Riemann surface of the model; different solutions lie on 
different sheets; there is room for solutions representing D-branes localized near different vacua 
by crossing the cuts that connect the sheets. 

To conclude the paper we would like to list a few open questions. The first concerns geometric 
transitions and gauge duals. Laufer's theorem, implies that only few smooth geometries can 
be interpreted as resolution of the singular conifold. These, as already observed in 0, correspond 
to asymptotically free gauge theories duals. Then, for the non Laufer's geometries one should 
formulate a definite gauge dual. Among these, one should find the non asymptotically free 
theories. It seems natural to guess the gauge dual of the geometries 0{—n) © 0{n — 2) to be 
then given hy M = 1 SYM with n + 1 Wess-Zumino multiplets in the adjoint representation of 
the U {N) gauge group and with superpotential given by the corresponding multi-matrix model 
one. Unfortunatly, we don't have convincing arguments to push further this hypothesis. 

Secondly, we elaborated a scheme which can be applied to Calabi-Yau manifolds of more 
general type than conifolds. It would be interesting to work out the effective coupling to the CY 
complex moduli of the reduced theory in such more general settings. 

The last comment concerns multi-matrix models. It is apparent that we have a method to 
solve any kind of two- or multi-matrix model with bilinear couplings. The next important step 
is to find analogous powerful tools to solve matrix models with more complicated couplings. 
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8 Appendix 

In this appendix we extend the method for the reduction of the holomorphic CS functional to 
the non compact CY geometry around a four cycle. This geometry, once one fixes the complex 
structure on the four manifold M describing the cycle, is fully determined to be total space of 
the canonical line bundle Km, that is the bundle of the top holomorphic forms on M. We denote 
this space as Xm = tot{KM)- 

Any atlas {U^a)} ^ extends to an atlas on Xm by = [/(q) x C. The complex manifold 
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is defined by the patching conditions 

Z(a) = f(a)(/3) (z(/3)) 

= [detX(^)(^)]-^P(^), where [X(„)(^)] = 5z(^)f(c,)(/3) (8.56) 

in any double patch intersection n In ()8.56|) and in the following, z = (z^ , z"^) denotes 
the two complex coordinates. The holomorphic (3, 0)-form on Xm is Q = dz^ A dz'^ A dp. 

Let us consider the topological B-model on Xm- In this case, D-branes can wrap the 4-cycle 
M and the theory describing the dynamics of these objects is obtained then by reducing the 
HCS functional to the D-brane world-volume. We consider here again only the case in which the 
gauge bundle E is trivial. 

The action of hCS is 

S(A) = — [ C, C = nATr(-AAdA + -AAAAA] (8.57) 
9s Jxm V2 3 y 

where A G r(O'i) (Xm). 

We split A = Azdz+Apdp and we set, because of the glueing prescriptions for the parallel and 
transverse components, A^dz = A^dz — ApT^pdz and Ap = Ap, where A = A^dz S T^^'^^ (M) 
is an anti-holomorphic 1-form on M, Ap G r{K^^) a section of the inverse anti-canonical line 
bundle and T^dz is the (0, 1) component of a reference connection on Km- 

The reduction prescription is that the matrix valued dynamical fields (A^jAp) are independent 
on the coordinate p along the fibre C. 

Direct calculation of the Lagrangian C in H8.57() for the above reduced configurations gives 

L = Cred = $1 A Tr Q 1^ A DAp + ApF^^^^'^ +AA TAp^^ dp (8.58) 

where D is the covariant derivative w.r.t. A. Notice that the above result does not depend on p. 

Introducing now a reference section K G r{KM Km), we define 0^^'°^ = KAp G r{KM) and 
fix the reference connection to be F = K~^dK. This way we get 

L = Cred = nK-' A Tr Q [a a D^^^''^ + </>(2,o)^(o,2)|^ ^- (8 

To reduce to a 4-form, we saturate the reduced Lagrangian with Kdp A dp so that the reduced 
HCS functional becomes just 

HCSred = -( Tr A 5<A(''0)+ 0(2,0)^(0,2)1 A 1 /■ ^ /^(2,0)^(0,2)\ (g gQ) 

9s Jm V2 U 9s Jm ^ ' 

which is the form used in |13j . 
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